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Abstract 

We use the procedure of pinched-weight Finite Energy Sum Rules (pFESR) 
to determine the OPE coefficients ae, • ■ ■ , cl\q of the flavor ud V-A correlator 
in terms of existing hadronic r decay data. We show by appropriate weight 
choices that the error on the dominant d = 6 contribution, which is known 
to be related to the K — > tttt matrix elements of the electroweak penguin 
operator in the chiral limit, may be reduced to below the ~ 15% level. The 
values we obtain for OPE coefficients with d > 8 are shown to naturally 
account for the discrepancies between our results for the d = 6 and d = 8 
terms and those of previous analyses, which were obtained neglecting d > 8 
contributions. 
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I. INTRODUCTION 



In a recent work [1], a pFESR analysis of the flavor ud two-point V-A current correlator, 
AI1(Q 2 ), was performed. This allowed extraction of the dimension six V-A OPE coefficient, 
a6, which is related by chiral symmetry to the K — > nn matrix element of the electroweak 
penguin operator, Q$. The result for led directly to an improved determination of e'/e in 
the chiral limit. 

The current paper is devoted to a more detailed account of this analysis. In it, we present 
the rationale for our choice of weight functions, describe the calculation of higher dimension 
OPE contributions, and discuss the relation of our work to previous treatments of the V-A 
correlator. Advantages of the particular version of the FESR formulation employed in our 
analysis are pointed out. 



A. Background 

We recall that non-strange hadronic r decay data provide access to the spectral functions 
of the flavor ud vector (V) and axial vector (A) current correlators [2-5]. With J^ A the 
standard V, A currents and the standard definitions of the spin J = 0, 1 parts of the 
correlators, 

i J d 4 xe^(0\T(.J^ A (x)J^ A (0y)\0) = (-(Tg 2 + gV) {<f) + <f<t ^ A (q 2 ) , (1) 
the ratios 

_ r[^~ -> z/ r hadrons v ,A (7)] , n x 

K V,A = ^7— — , . , (/J 

(with (7) indicating additional photons or lepton pairs) are expressible as weighted integrals 
over the corresponding spectral functions py\ = - Imlll/^. Working with the combinations 

IlyJ^(s) = II^(s) + II^(s) and sIly^ A (s), which have no kinematic singularities, one has 
explicitly [2-5], 



Rv,a = 12tt 2 S E w \V u d\ 2 f — 5- (l - 

Jo mi \ 



t ds ( _ s 



(3) 



where Sew = 1.0194 ± 0.0040 represents the leading electroweak corrections [6] 1 , and V ud 
is the ud CKM matrix element 2 . Since the integrals over the J = part of the spectral 



X A recent update [7] yields Sew = 1-0201 ± 0.0003, compatible, within errors, with the value 
1.0194 ±0.0040 quoted above, and employed in our recent paper [1]. Since the ±0.0040 uncertainty 
on Sew produces a negligible contribution to our total errors below, we have chosen to retain the 
input employed in Ref. [1] in what follows. 

2 The additional radiative correction, conventionally denoted S' EW , has been dropped in writing 
this equation since it cancels in the V — A difference which is the subject of this paper. 
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function are saturated by the pion pole contribution, up to numerically negligible corrections 
of 0(m^ d ), non-strange hadronic r decay data provide detailed information on the sum 



(o 
Pv 



s) + Pa +1> (s). For states containing only pions, G-parity allows an unambiguous 
separation of the V and A components of this sum. In the range where the decay to states 
containing kaon pairs is negligible (say s < 2 GeV 2 ), the individual V and A terms, and hence 



also the difference Ap = p y v 



.(o+i) 
Pa > 



are thus known very accurately from experiment. 



Knowledge of the V and A spectral functions allows access to the corresponding corre- 
lators through the use of either dispersion relations or FESR's. The latter may be taken to 
have the form 



/ dsp(s)w(s) —- — : & dsH(s)w(s) , 

J 8 t h 27T2 J|s|=s 



(4) 



valid for any w(s) analytic in the region of the contour, and any H(s) without kinematic 
singularities. An example is the standard OPE representation of Ry,A [2-5], 



Rv,a = QttSew \V ud \ 2 i f 

J\s\-- 



ds 

%l ml 



mi 



1 + 2- 



mi 



m\ ' 



, (5) 



which results from the application of Eq. (4) to Eq. (3). 

If one works at sufficiently large s that the OPE representation of H(s) may be used 
reliably on the RHS of Eq. (4), appropriate choices for the weight w(s) allow one to deter- 
mine OPE contributions of different dimension, d, in terms of experimental data for p(s) 
(analogous statements are true for the corresponding dispersion relations and/or their Borel 
transforms). To reflect the fact that II(s) will differ from its OPE representation over at 
least some portion of the contour \s\ = sq, we recast Eq. (4) in the general form 



/ dsp(s)w(s) + R[s ,w] = — — I dsU OPE (s)w(s) 

Js th l-Kl J\s\=sq 



where 



R[s ,w] = —A / ds (IIope(s) - n(s)) w(s) . 

AlXl J\s\=s 



(6) 



(7) 



R[so,w] then quantifies what is usually referred to as "OPE breakdown" or "duality viola- 
tion". 

Due to the intractibility of strong-coupling QCD, it is not possible at present to obtain 
an analytic expression for R[so, w], and its neglect (common to all FESR analyses) therefore 
represents a key dynamical assumption. There exist several strategies, however, to minimize 
the impact of this assumption: 

• Work at the highest possible s = s max , where s max is the maximum value of s for 
which p(s) is experimentally known. At high sq, one has increased confidence in 
the reliability of the OPE. One can check the stability of any nominal OPE output 
against changes in s to assess theoretical uncertainties associated with possible OPE 
breakdown. 
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• Work in the vicinity of certain 'optimal' s values s , called duality points [9]. The set 
of all such Sg^ is nothing but the zeros of R[sq,w] for certain special weight choices. 
By "special" is meant that the zeros can be determined independent of the values of 
any unknown OPE condensates. For example, the weights w(s) = 1 and w(s) = s are 
special for the V-A correlator, since the d — 2 and d — 4 OPE contributions are known 
to be zero in the chiral limit (this is the OPE statement of the two Weinberg sum rules 
(WSR's) [10]). In general, the sets of zeros of R[sq,w\ for different correlators and/or 
different weights are different 3 . The duality point approach relies on the observation 
that, for the V-A correlator, the w(s) = 1 duality points lie close to the corresponding 
w(s) = s duality points (two such points exist in the interval < s < m%.). This 
is taken to suggest the possibility that the zeros of R[so,w] for all w(s) might be 
(approximately) the same. If so, sum rules based on other weights would be reasonably 
satisfied at the w(s) = 1, w(s) = s duality points. Such sum rules, restricted to these 
values of Sq, could then be used to extract unknown OPE coefficients. However, the 
uncertainty about how close the true duality point for a given sum rule is to that for the 
w(s) — 1, s sum rules will produce a corresponding uncertainty in the extracted OPE 
coefficients. This uncertainty can be large if the sq dependence of the corresponding 
spectral integrals is strong (as it is, for example, for the weights w(s) = s k with k > 2). 

• Work with 'pinched weights', i.e., those satisfying w(s ) = 0. Such weights suppress 
OPE contributions from the region of the contour near the timelike real axis where 
[Il0P£ — II] is expected to be largest [12]. R[sq,w] will then be small at those scales, 
so, for which the region of OPE breakdown, on the contour \s\ = s , is restricted to the 
vicinity of the timelike point. If the region of scales for which this is true extends down 
as far as the experimentally accessible region, one can find a window of s values within 
which the data-based spectral integrals admit an OPE-like representation (as given by 
the RHS of Eq. (6)) for suitable choices of the unknown QCD parameters (i.e. the 
appropriate OPE condensates). A successful OPE/data match implies that R[s ,w] is 
not detectable, within experimental errors, in the given analysis window. The analysis 
window then represents an extended duality interval, since every point in it, in the 
sense of the terminology above, is a duality point. Moreover, working with weights 
which differ significantly in the ways they weight the experimental spectral data, but 
whose integrated OPE contributions involve the same set of QCD parameters, allows 
one to perform additional checks. 

We shall follow the last of these approaches and employ pFESR's to analyze the ud V-A 
correlator. Evidence in support of this choice can be inferred from the results of FESR 
studies of the flavor ud V and A correlators, where one knows with good accuracy both 
the data and the OPE integrals for s above ~ 2 GeV 2 . These studies show that FESR's 
based on the unpinched weights, w(s) = s k (k = 0, • • • , 3), are rather poorly satisfied over 
the range 2 GeV 2 < sq < m 2 T [11,16] (i.e., at these scales, R[so,s k ] is typically large). In 
contrast, the FESR predictions, Eq. (5), for R v and Ra, which are obtained by taking 



3 See, for example, Figure 1 of Ref. [16]). 
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the appropriate linear combinations of the w(s) = l,s 2 and s 3 FESR's, with sq set equal 
to m 2 . and R[m 2 ,w] to zero, are in extremely good agreement with experiment [18,19]. 
The failure of the s fc -weighted FESR's is a manifestation of the breakdown of the OPE 
representation near the timelike real s axis for insufficiently large s , as shown by Poggio, 
Quinn and Weinberg (PQW) [12]. The success of the OPE predictions for R v ,a presumably 
arises from the suppression of this danger region by the (double) zeros of the kinematic 
weights at s — m 2 (the edge of hadronic phase space). It turns out that for any weight of 
the form either u>jv(y) = (1 — y)(l + Ay) or wd(u) = (1 — 2/) 2 (l + Ay) (with A arbitrary 
and y = s/so), the corresponding i?y,A-like pFESR is extremely well satisfied for all so 
in the range 2 GeV 2 < s < m 2 . [11,16]. This indicates that for the separate V and A 
correlators, and for such "intermediate" scales, the OPE breakdown is closely localized to 
the vicinity of the timelike real axis. At these scales, it appears safe to neglect R[sq,w] 
also for other correlators provided w(y) satisfies w(y = 1) = 0, but not otherwise. A more 
detailed discussion of these issues may be found in Section V. 



B. Summary of Content 

In this paper we focus on the difference An=IIy )+1 ' ) — U^ +1 ^ of the J = + 1 components 
of the flavor ud V and A correlators. In the chiral limit, its OPE is purely non-perturbative, 
with contributions beginning at dimension d — 6. The smallness of m U ( i means that the 
physical OPE will be dominated by d = 6 (and higher) terms, at least until one gets to 
extremely large scales. Accurate data for the associated spectral function Ap thus allow 
the extraction of various vacuum condensate combinations. Of particular interest are the 
two condensates appearing in the d = 6 part of [An] OPE , which turn out to determine the 
chiral limit values of the K — > tttt matrix elements of the electroweak penguin operators 
Q7,8 [13]. The d > 6 terms in [An] OPE , which enter dispersive sum rules for these matrix 
elements [14,15], are also of phenomenological interest since a determination of their values 
would allow the dispersive determination of Ref. [15] to be performed at lower scales, where 
uncertainties associated with the classical chiral sum constraints are drastically reduced. 

We will extract the higher dimension (d > 4) terms appearing in the OPE of All by con- 
structing a set of pFESR's designed in such a way as to minimize the impact of experimental 
errors. As we will show below, it is possible to make such determinations for d — 6, • • • , 16 
with good accuracy using the existing experimental r decay data base. 

In Section II, we detail the input required for the OPE and data sides of the flavor ud V-A 
r decay sum rules and discuss some practical considerations relevant to the choice of pFESR 
weight. In Section III, we describe how, by appropriate pFESR choices, it is possible to (1) 
significantly improve on previous determinations of the d — 6 and d = 8 OPE contributions 
and (2) at the same time, extract OPE contributions with dimensions d — 10, • • • , 16 not 
obtained in those earlier analyses. A comparison with previous analyses is presented in 
Section IV. An expanded discussion of the issue of duality violation is given in Section V. 
This section also contains an outline of the techniques we have employed to test for the 
presence of possible residual duality violation in our analysis of the V-A correlator. Certain 
details of these tests which are relevant to the comparison to earlier work are deferred to 
the Appendix. Our conclusions, together with a brief discussion, are given in Section VI. 
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The implications of our results for the chiral-limit values of the electroweak penguin matrix 
elements have already been worked out in Ref. [1]. 

II. SUM RULE ANALYSES OF THE V-A CORRELATOR An 

We now describe the input required for the OPE and data sides of the flavor ud V-A sum 
rules. The emphasis is on practical considerations relevant to the choice of pFESR weight 
functions. 

A. The data side 

We shall employ both the ALEPH and OPAL data for Ap [18,19]. The respective ALEPH 
and OPAL spectral functions are displayed in Fig. 1 4 . 

In the case of ALEPH, we use the publicly available data files corresponding to the 
1998 analysis, whose overall normalization was set by the preliminary result for the rescaled 
total strange hadronic branching fraction, R us = B us /B e , R us = 0.155, and the 1998 PDG 
values of B e and B^. The 1999 published version, R us = 0.161 [20] and the recent update, 
R us = 0.1625 [21,22] both differ slightly from the preliminary value. This change, together 
with minor changes in the values of B e , B^ and the r — > itv T branching fraction, B n , 
necessitates a small global rescaling of the 1998 V-A data and covariance matrix 5 . An input 
value of V U( i is required to convert from the experimental number distribution provided by 
ALEPH to Ap. We have taken this to be V u d = 0.9742 ± 0.0016, a value which spans both 
the range based on the K e3 decay analysis and that based on the combination of + — > + 
nuclear decays and neutron decay, as quoted in PDG2002 [23]. In the case of OPAL, we 
use the publicly available data files for Ap and its correlation matrix, corresponding to the 
results of Ref. [19]. These files were constructed using a central value V u d = 0.9753. We 
have, therefore, performed a small global rescaling in order to work with ALEPH and OPAL 
versions of the spectral function which both correspond to V ud = 0.9742 ± 0.0016. 

As noted previously, the data are very accurate below s ~ 2 GeV 2 . Near the kinematic 
endpoint y T = 1 (y T = s/m%), however, the errors on Ap become large. This is a consequence 
of several factors: 

1. The event rate becomes small in that region due to phase space suppression. 

2. There is, at present, no complete separation of V and A contributions to the spectrum 
for states containing a KK pair and > 1 7r's. 



4 The points shown represent bin-averaged Ap(s) values and are plotted at the midpoints of the 
ALEPH/OPAL experimental bins. 

5 We thank Shaomin Chen for pointing out the necessity of this rescaling to us. For R us = 0.1625, 
the PDG2002 average value [23] B e = 0.1781, as implied by p-T universality and the 7r M2 value 
of F n , the rescaling turns out to be 1.003, i.e., very close to 1. 
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3. In order to extract Ap from the experimental decay distribution (c./. Eq. (3)) one 
must divide by the kinematic weight factor w T (y T ) = (1 — y T ) 2 {l + 2y r ). The double 
zero of w T at s = m 2 thus amplifies the errors on the V-A number distributions for 
those s near s = m 2 . 

In view of item 3, pinched weights with only a single zero at s = Sq will weight the ex- 
perimental number distribution and errors with a factor which diverges as s — > s for s 
near m 2 . Such behavior is to be avoided if one wishes to keep the errors on the weighted 
spectral integrals under control. For this reason we restrict our attention in the following 
to pFESR's based on polynomial weights of the form p(y)(l — y) 2 (y = s/sq). Though this 
restriction is forced on us by necessity, it has the virtue of enforcing a stronger suppression 
of OPE contributions from the vicinity of the timelike real axis, and hence of improving the 
reliability of the OPE side of the pFESR's. 

An important practical consideration in choosing pFESR weights is the non-positive 
definiteness of Ap. Even with the very precise data below 2 GeV 2 , weighted integrals 
which involve significant cancellations between contributions from the regions of positive 
and negative Ap will have much larger fractional errors than would be expected based only 
on the accuracy of the spectral data alone. For some of the weights employed previously 
in the literature, for example, the V-A cancellation is at the level of a few percent of the 
individual V and A integrals, leading to large errors and significant sensitivities (as large as 
20%) to the exact treatment of the n pole contribution. Avoiding strong cancellations of 
this type is crucial to reducing the errors on the final determinations of the various vacuum 
condensates. To quantify this point in our discussions below, we introduce a quantity ry-A 
defined as the ratio of the V-A spectral integral to the corresponding vector spectral integral. 



B. The OPE side 



The OPE representation of All is schematically of the form 

AIT(Q 2 )= E§- (8) 

d=2,4,- V 

where Q 2 = — s. Perturbative corrections lead to logarithmic dependences of the X d on Q 2 . 
To NLO in QCD one has: 

X d = a d (fi) + b d log (^Pj . (9) 

The bd(fi) are known explicitly for d — 2, 4, 6, but not for higher d. For polynomial weights, 
OPE contributions proportional to bd(fj) involve the integrals Jj x i =1 dx x k log(x). For the 
weights w(y) = p(y)(l—y) 2 employed in our analysis the variations of sign in the coefficients 
of w(y) produce significant cancellations (and hence additional numerical suppressions) of 
these contributions relative to those of the leading non-logarithmic terms. We thus consider 
it very safe to follow earlier analyses in neglecting such corrections for d > 8. The remaining 
non- logarithmic OPE integral contributions follow from 
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-1 

2ni J\s\=so 



ds 



Q d 



h,(d/2)- 



Od 
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(10) 



Weights of degree 2 thus contain leading OPE contributions up to d = 6, those of degree 
3 contributions up to d = 8, etc.. Neglect of d — 10, • • • , 2N + 2 contributions in pFESR's 
based on weights with degree N > 3 is therefore dangerous unless one is working at Sq large 
enough that such contributions may be taken to be safely small. Typically one does not 
know a priori how large an s is "large enough"; however, the stronger l/s -dependence of 
the higher d integrals allows this question to be addressed post facto, provided one works 
with a range of s large enough to expose the presence of higher d contributions which 
may have been omitted when they should not have been. If one finds that the range of So 
employed is such that the presence of such contributions is indicated, one can use the pFESR 
in question to place constraints on the relevant higher dimension ad terms. Obviously, both 
the reliability of the post facto check and the accuracy of the higher d extraction will be 
enhanced for pFESR's having fewer separate a d contributions on the OPE side of the sum 
rule and larger separations between the dimensions of the contributions which do occur. 

The d = 2 term in Eq. (8) is of 0(m 2 u d ). That it can be safely neglected can be confirmed 
numerically by integrating the J = + 1 expression of Ref. [25], which is known to 0(a 2 s ). 6 

The d = 4 term in Eq. (8) is given by [3,29] 



[An(Q 2 )] 



d=A 



8_ 59_ 2 \ ((m u + m d )uu) 



i a + Y a ~ 



Q 4 



(ii) 



where a = a s (Q 2 )/ir, with a s (Q 2 ) the running coupling at scale fj 2 = Q 2 in the MS scheme. 
The quark condensate factor can be evaluated using the GMOR relation [30] 



((m u + m d )uu) = -F 2 ml 



;i2) 



which is accurate to better than 6% [31]. We compute the weighted integrals of [AIl(Q 2 )] d=4 
using the "contour improvement" scheme [4,32], taking for a the version corresponding to 



6 The reader might worry that the rather bad behavior of the integrated J = 0, d = 2 OPE 
series precludes reliably subtracting the non-7r-pole part of the J = contribution from the data, 
and hence prevents us from making such a definitive statement. While it is true that (i) for 
the kinematic weight case shown above, the last three terms in the integrated J = d = 2 series 
(which is known to 0(a%) [26]) are actually increasing [26,27], even at the scale sq = m 2 and (ii) the 
0(«s) -truncated J = d = 2 OPE integrals corresponding to different u (k,0) spectral weights" [4] 
display a significant unphysical dependence on k [28], this turns out not to be a problem. The 
reason is that the behavior of the integrated J = series has been investigated in the analogous 
case involving the flavor us currents, where additional sum rule constraints were shown to allow 
a determination of the corresponding spectral integrals [28]. The 0(a\ )-truncated J = OPE 
estimates were found to represent significant overestimates [28]. We may thus use the (albeit 
poorly behaved) OPE determinations to conclude that, apart from the tt pole contribution, the 
d = 2 J = contributions to the measured spectral distribution are indeed completely negligible. 
The J = + 1 part of the spectral function can thus be reliably determined. The integrated series 
in a s for the d = 2 J = + 1 OPE contribution converges well, and is numerically negligible. 
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4-loop running [33] with the ALEPH determination [18] a s (m 2 T ) = 0.334 ± 0.022 as input. 
This contribution represents only a small correction to the dominant d = 6 term because 
of the 0{m Ut< j) chiral suppression. In the numerical analysis we have expanded to ±20% 
the errors assigned to the GMOR evaluation of the d = 4 OPE contributions in order to 
account for the truncation of the series for the Wilson coefficient at O (a 2 ). Because the 
d — 4 contribution is so small, the resulting contribution to the total error is, however, 
negligible. 

Observe that in the chiral limit the d = 2 and d = 4 contributions are zero. Taking 
0,2 = CL4 = is then the OPE implementation of the first and second WSR's. To the extent 
that we use a 2 = and a 4 ~ m 2 the WSR's are built into our procedure. 

For the d = 6 contribution, there exist several determinations in the literature [15,34,35], 
corresponding to different schemes for the choice of evanescent operator basis [36] . Since one 
of our goals is to use our results for the d = 6 contribution to improve the determination of the 
chiral limit value of the electroweak penguin contribution to the K — > Tin decay amplitudes, 
we employ the most recent determination [15], which corresponds to the same scheme as 
used in the calculation of the Wilson coefficients of the effective weak Hamiltonian [37] 7 . 
To simplify the later application of our results it is also convenient to work with the vacuum 
condensates (Oi) and (0 8 ) defined in Ref. [15], 



^^q - <?7^75yg qTi^q) 



7-3 

(Oi) = {qj^—q <rt 
(Os) = (qi^q ll^q ~ qi.l^q qi'l^q) 



(13) 



where q — u,d, s, r 3 is a Pauli (flavor) matrix, and {A a } are the Gell Mann color matrices. 
With these choices one has 



Q 2 

cfe(A0 + ^e(A*)ln — 



(14) 



with 



ae(pL) = 2 [27r(a s 8 ) M + A s (a 2 s 8 )^ + Ai(a^Oi)^ 
be(fi) = 2 \B 8 (a 2 s 8 )^ + B 1 (a 2 1 )^ 



(15) 



where A ± , A 8 , B 1 and B 8 are the coefficients tabulated in Ref. [15] 8 . They depend on the 
number of active flavors, the scheme employed for 75, and the evanescent operator basis. 
For Nf = 3, the values for the NDR and HV 75 schemes are 



7 An independent determination of [AIl] rf=6 in this scheme was given in Ref. [38]. The results 
quoted in version 2 of this reference are now in agreement with those of Ref. [15]. 

8 The coefficients and be appearing here differ from those of Ref. [15] by a factor of 2. This reflects 
the fact that in Ref. [15] the coefficients correspond to the neutral isovector current correlator, while 
here they correspond to the charged isovector current correlator. The (isospin) factor of 2 has been 
made explicit in Eqs. (15). The A\, A 8 , B\ and B 8 of Eqs. (15) thus have the same numerical 
values as in Ref. [15]. 
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A ± = 2 (NDR), -10/3 (HV) 
A 8 = 25/4 (NDR), 21/4 (HV) 
B 1 = 8/3 (NDR and HV) 

B 8 = -1 (NDR and HV) . (16) 

The logarithmic (Bi >8 ) terms turn out to play a very small role in the analysis, though 
we have kept them for completeness. We do this by first writing b$ = a^ib^/ao) and then 
employing the existing dispersive determination of (Oi) and (Og) [15] to estimate r 6 = b^/a^. 
With this estimate as input, the integrated d = 6 OPE contribution is now, like the non- 
logarithmic d = 6 term, proportional to a 6 . The overall a 6 factor multiplying the full d = 6 
contribution is then to be fit to data. The central value for r 6 turns out to be very small, 
~ —0.03. Since only the first term in the expansion of b$ in powers of a s is known, we assign 
a (conservative) 50% uncertainty to this estimate. 
For d = 8 and higher we take 

{An(Q% = ^. (17) 

With this notation, a 8 is identical to (0 8 ) of Ref. [8] and 8 of Ref. [24]. It is also twice 
the negative of the integral M 3 of Ref. [38,39], independent of s , in the absence of duality 
violation in the s 3 -weighted FESR. 



III. EXTRACTION OF OPE CONDENSATES FROM PFESR'S 
A. Choice of pFESR Weights and the so Analysis Window 

We consider a sequence of pFESR's designed to simplify the extraction of the OPE 
coefficients of An. Working with pFESR's allows us to take advantage of the freedom in 
the choice of weight profile, and hence, by construction, to avoid strong V-A cancellations. 
The freedom of weight choice also allows us to considerably simplify the task of constraining 
the d > 8 contributions and separating them from the d = 6 and d = 8 contributions. 

Within the space of pFESR weights w(y) = p(y)(l — y) 2 employed in our analysis, the 
weight of lowest possible degree is w(y) = (1 — y) 2 . In a zero-error world, the corresponding 
pFESR would allow an extraction of a 6 . Unfortunately, this weight produces a high degree of 
V-A cancellation, and hence is not practical for use when employed with present experimental 
data. We thus consider weights of degree 3 (the highest degree possible involving no a d 
contributions with d > 8), w(y) = (1 — y) 2 (l + Ay). There will be some value of A for which 
the fractional errors on the spectral integrals are minimized. It turns out that this value is 
almost exactly equal to —3. The pFESR based on 

Wl (y) = (l-yf(l-3y) (18) 

will then provide the most restrictive constraint on a^, a§, and this is our first choice of 
weight. The weight w 2 (y) = y(l — y) 2 also has reduced V-A cancellation, and provides 
independent constraints on the determination of a 6 and a 8 , since degree three polynomials 
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yield only 04, a$ and a§ OPE contributions, and 04 is small. The so dependence of both 
the W\- and u^-weighted spectral integrals will be well-described using only two parameters, 
a 6 and a 8 , provided the use of the OPE representation is justified. w 2 has been chosen, by 
construction, to weight Ap(s) = Ap(s y) with a profile very different from u>i to make this 
test of the reliability of the OPE representation as non-trivial as possible. 

We determine our so analysis window by fixing the upper edge at sq = 3.15 GeV 2 ~ m 2 
and decreasing the lower edge until the fitted coefficients cease to be consistent (within 
experimental errors). Since a 6 is the most accurately determined coefficient we use it as our 
basic monitor of the onset of duality violation. We find that duality violation for the V-A 
correlator and the Wi, w 2 weight set begins to set in below s ~ 1.8 GeV 2 , and hence we fix 
the lower edge of our analysis window at 1.95 GeV 2 . 

To investigate d > 8 contributions it is convenient to construct weights for which the 
only OPE contributions with d > 4 are those proportional to a 6 and aa, with d = 10, 12, • • •. 
The possibility of working with s down to ~ 2 GeV 2 is also helpful since an increased 
range of s creates an increased variation in the relative size of the d = 10, 12, • • • and 
d — 6 contributions over the analysis window, and hence improves our ability to perform the 
separation of contributions of different dimension. Weights having a double zero at y — 1, 
reduced V-A cancellation, and only a single contribution beyond d — 6, are 



wn(v) = y 



1 - { W^i) y+ {w=^ ]y 



N 



N = 2,3,4,5,6 . (19) 



The overall factor of y has been introduced in order to reduce the level of V-A cancellation. 
The case N = 2 corresponds to the previously introduced weight w 2 = y{l — y) 2 . For N > 2, 
w N(y) produces contributions proportional to 04, a§ and a 2 N+i on the OPE side of the sum 
rule. We consider N up to 6, and hence ad contributions with d up to 16. 9 Since each 
of the resulting sum rules allows a determination of both a§ and a 2 N+4, the consistency 
of the a,Q solutions obtained from the W\ through wq pFESR's also provides a strong self- 
consistency constraint on the reliability of the analysis. Further constraints on a§ and the 
higher dimension can be obtained by considering the weights 

w, +N (y) = y [1 - (^2) y 2 + (j^) v N ] N = 3 > 4 > 5 > 6 ■ ( 20 ) 

These weights produce a 4 , a 8 and ad contributions with d = 10, 12, 14 and 16 for A^ = 3, 4, 5 
and 6, respectively. The as and ad>8 values extracted using wj through wio should be 
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N > 2, 



We also investigated pFESR's based on the weights wjsr(y) = 1 — ( jy-i ) D (iv^T 
which produce only and a 2 N+2 OPE contributions. The N = 2 case is just w\{y) = (1 — y) 2 . 
For larger N the smallness of the d = 4 contributions would, in principle, make pFESR's based on 
these weights good choices for determining the higher dimension terms. The V-A cancellations 
for the Wn family are, however, considerably stronger than for the Wn family of Eq. (19), making 
the errors on the extracted significantly larger than those obtained using the pFESR's based 
on w\ through w\q. While the results for the obtained using the two sets of sum rules are in 
excellent agreement, the larger errors make the analysis based on the wn inferior to that based on 
the wn, at least with current experimental data as input. 
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consistent with those obtained using w\ through w^, provided the OPE representation of 
All is reliable for the Sq employed in our analysis. We find the consistency is excellent for 
all the a d with d > 4. 

Finally, we observe that pFESRs based on the weights of Eqs. (19) and (20) allow one in 
principle to extract condensates of even higher dimension. With the present experimental 
errors, however, higher degree pFESR's effectively work with a smaller analysis window, 
localized around sq = 2 GeV 2 (points at higher sq suffer from much larger experimental 
errors, and become irrelevant in the analysis). This feature weakens the power of this 
method to detect inconsistencies through the use of an extended s analysis window. We 
therefore quote our results for the condensates only up to d — 16. 

With the above choice of weights, and assuming R[so,w] = 0, the w\ through wio 
pFESR's may be written as 

Jw n (s ) = f Wn ({ a d}; s ) 



where 



J Wn (s ) = f dsw n (— ) Ap(s) + — I dsw n (—") [AIIope 

JO \Sq/ 2m J\s\=s VSn/ 



fu, n ({ a d};s ) = 



s 

ds w. 



2m J\s\ 

The explicit form for the OPE integrals is: 



|s|=so 

[An PE(s)] (i>4 



»] 



d=4 



f Wl ({ad}] so) = -2 a 6 



fw 2 {{0-d}] So) 
fw 3 ({ a d}', So) 
fw 4 (i a d}; So) 

f W5 ({ad}; so) 
f W6 {{a d }; s ) 
f W7 ({ad}; s ) 
f W8 ({a d }; s ) 
fw 9 ({a>d}; so) 
fw w ({ad}; so) 



s 2 o 



_6 

Wo 

3 



1 + r 6 log y 
1 + r 6 log | 
1 + r 6 log 
1 + r 6 log 
1 + r 6 log 



So 

/A 



3 

+ iT 6 



+ 



3a 8 



So 



£o 



J 1 



£o 



2 a 6 



1 + r 6 log 



'sq S 



o^Clo + 



3<2fi 



+ 



s Sq 

8 2a 8 

— 2 r G a 6 + — 3~ 

s 
5a 8 



2aio 

s o 
ai2 



s 

1 

+ 2 r6 



+ 3 r6 
3 

4 

+ ~r 6 

5 



aio 
24 

Ql2 

3s 5 



4s| 

Ql6 

5sl 



3ag 
5 

~2?o r6ae + 34 + 
12 3a 8 



„5 
s 

2a 



14 



3sg 
ai6 



5Sq 2Sn 2Sn 



(21) 

(22) 
(23) 



(24) 



'o -"o 

Note that the small, known d — 4 OPE contribution has been moved to the spectral integral 
side in defining J TO „(so). 
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B. pFESR Fit: Input and Results 



For our final results, we proceed in two steps. In the first step, we extract a preferred 
value for a 6 in an analysis employing only the weights W\ and w 2 [1]. Such an analysis 
is "maximally safe" in the sense that the numerical suppressions of the integrated d > 8 
OPE logarithmic corrections are strongest when d — 2n, where n is the degree of the pFESR 
polynomial, is as large as possible; the neglect of such d > 8 logarithmic terms in the OPE 
is thus safest when one uses the weight (s) of the minimum possible degree. As explained 
above, the accuracy of current data means that the lowest such degree which still allows 
an accurate extraction of a 6 is 3. In the second step, we perform a combined least-squares 
fit for the coefficients a^, - ■ ■ , using, for each of the weights w± through tuio, defined in 
Eqs. (18), (19), and (20), the set of 7 s values 1.95 + 0.2A; GeV 2 , k = 0, • • • , 6, which span 
the range from s ~ 2 GeV 2 to 3.15 GeV 2 ~ m 2 . 

On the data side we use as input for the analyses based on both the ALEPH and OPAL 
data B e = 0.1781 ± 0.0006 [23], F w = 92.4 ± 0.07 ± 0.25 MeV [23], S EW = 1.0194 ± 0.0040, 
and \V ud \ = 0.9742 ± 0.0016. The rescaling of the 1998 ALEPH data is determined using 
Rud = B ud / B e = 3.480 ± 0.014 [22]. This value is based on the most recent update, R us = 
0.1625 ± 0.0066 [21,22], in combination with the PDG2002 average for B e (quoted above), 
and the assumption of /i-e universality. On the OPE side we use ((m u + m d )uu) = —F%m^, 
and r 6 = b 6 /a 6 = -0.030 ± 0.015. 

In listing final errors for the ALEPH-based analysis we quote separately the errors pro- 
duced by the uncertainties in the ALEPH number distribution, and those due to all other 
sources, including the uncertainties on the OPE input quantities and r 6 . The former 
are calculated using the rescaled ALEPH covariance matrix. The latter are combined in 
quadrature. 

In the analysis based on the OPAL data, we again quote two uncertainties. The first is 
that computed using the OPAL covariance matrix, the second that obtained by combining 
in quadrature the errors associated with uncertainties in all other input parameters (V u d, 
Sew, o 4 , and r 6 ). 

Fits to the ALEPH data 

The results of the "maximally safe" analysis for a 6 and a s are 10 



10 Due to strong correlations between the data integrals for different sq and different weights, the 
fit values are obtained by minimizing the sum of the squared deviations between the data and OPE 
integrals, weighted by the inverse of the diagonal elements of the covariance matrix for the set of 
data integrals [1]. With this procedure, as is well known, the one-sigma errors and rms errors do 
not coincide. The former are smaller, and underestimate the variation in the fitted a d produced 
by variations in the input experimental data. All errors quoted in what follows are, therefore, the 
(larger) rms errors, i.e., the square roots of the diagonal elements of the covariance matrix for 
the {ad} solution set. The fitted values are, of course, also strongly correlated, and it is crucial 
to employ the full covariance matrix for the solution set if one wishes to have accurate errors for 
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a 6 = - (4.45 ± 0.61 ± 0.34) x 10~ 3 GeV 6 

a 8 = - (6.16 ± 2.78 ± 1.40) x 10~ 3 GeV 8 . (25) 
For the "combined fit" analysis, we find 



a 6 = 


-(4.54 ±0.83 ±0.18) x 10~ 3 


GeV 6 


a 8 = 


- (5.70 ±3.72 ±0.64) x 10" 3 


GeV 8 


a w = 


(4.82 ± 1.02 ±0.20) x 10~ 2 


GeV 10 


O.Y2 = 


-(1.60 ±0.26 ±0.05) x lO^ 1 


GeV 12 


ai4 = 


(4.26 ±0.62 ±0.14) x 10" 1 


GeV 14 


a i6 = 


-(1.03 ±0.14 ±0.03) GeV 16 





(26) 



Fits to the OPAL data 

The results of the "maximally safe" analysis for a e and a 8 are 



a 6 = - (5.43 ± 0.72 ± 0.25) x 10" 3 GeV 6 

a 8 = - (1.35 ± 3.32 ± 1.0) x 10~ 3 GeV 8 . (27) 
For the combined analysis, we find 

a 6 = - (5.06 ± 0.89 ± 0.12) x 10 -3 GeV 6 

a 8 = - (3.12 ± 3.82 ± 0.45) x 10 -3 GeV 8 
a w = (3.87 ± 1.06 ±0.10) x 10~ 2 GeV 10 
a 12 = - (1.32 ± 0.27 ± 0.03) x 10" 1 GeV 12 
a M = (3.54 ± 0.66 ± 0.06) x 10" 1 GeV 14 

a 16 = _ (0.85 ± 0.15 ± 0.02) GeV 16 . (28) 

We note that the ALEPH and OPAL determinations of OPE coefficients are in good 
agreement within errors. There is also extremely good agreement between the combined-fit 
and maximally- safe-fit values for a 6 and a 8 in both the ALEPH and OPAL cases, providing 
strong post facto support for the neglect of the higher d logarithmic corrections. One further 
point of relevance to the self-consistency of the analysis, not evident from the results quoted 
above, is the following. For each of the ten pFESR's considered above, it is possible, because 
of the different so-dependence of contributions of different dimension, to extract values for 
the two unknown (d > 4) a d coefficients occuring on the OPE side of the sum rule in question. 
One can then compare the values of a given obtained using various different individual 
pFESR's. It turns out that the agreement among the results of different single-pFESR 



various sums of higher d OPE contributions such as those that enter the dispersive test of the 
solution set described in the Appendix, or those required if one wishes to perform the residual 
weight analysis for the K — ► tttt EW penguin matrix elements at lower scales [15]. 
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analyses is excellent for all the a d , d — 6, ■ ■ ■ , 16. By construction six such determinations, 
and hence six such consistency tests, exist for a§ and a 8 . 

Our combined fit leads to a determination of the six parameters a 6 . . .a 16 . Because of 
the strong correlations between data integrals corresponding to different s and/or different 
weights, the resulting fit parameters are highly correlated. If Cdd 1 is the DD' element of 
the correlation matrix, we find that the smallest of the |Cdd'| is 0.90 for the solution set 
associated with the ALEPH data and 0.82 for that associated with the OPAL data. The 
full covariance matrices are available upon request. 

C. The Optimized OPE/Spectral Integral Match 

It is important to verify that, after fitting the OPE coefficients a d) the resulting OPE 
integrals f w „({a d }; s ) provide a good match to the corresponding spectral integrals J Wn (s ) 
over the whole of the Sq analysis window. Failure to achieve such a match would represent a 
clear sign of duality violation. In Figs. 2, 3 and 4 we display the quality of the f Wn {{dd}', so) 
/ Jw n {so) match for the combined fit to the ALEPH data. (The match for the combined fit 
to the OPAL data is of identical quality, and hence not shown separately.) Fig. 2 shows the 
results for the w 1 and w 2 pFESR's 11 , Fig. 3 for the w 3 through w 6 pFESR's, and Fig. 4 for 
the w-i through u> 10 pFESR's. Our results for f Wn {{ a d\\ So)> corresponding to Eqs. (26), are 
given by the solid lines. There is clearly no sign of duality violation for any of the pFESR's 
employed at any of the scales, sq, in our analysis window. Improved data would reduce the 
errors on J Wn (s ) and allow us to sharpen this test even further. Also shown for comparison 
in each figure are the OPE results corresponding to the a 6 , a 8 fits of Refs. [8,24], where, 
as in those references we take as central input values a d = for d > 8. The inclusion of 
d > 8 contributions clearly leads to a significantly improved fit to the data, as well as a 
significantly reduced error on the determination, in particular, of a§. 

The excellent agreement between the optimized OPE representation and the correspond- 
ing data integrals displayed in Figs. 2 through 4, while a necessary condition that significant 
duality violation be absent from our analysis, is not a sufficient one. In order to investi- 
gate this question further, we have performed a number of additional tests on our solution 
sets. Since several of these tests correspond to sum rules studied in earlier analyses of the 
V-A correlator, we first discuss the relation between our results and those of these earlier 
analyses. Having introduced the relevant sum rules as part of this discussion, we will then 
return to a discussion of the additional tests which such sum rules allow us to perform on 
our solution sets in Section V. 



We plot only the results of the combined fit in this case since they are indistinguishable from 
those of the "maximally-safe" fit on the scale of the figure. 
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IV. PREVIOUS ANALYSES 



Several determinations of the d — 6 and d = 8 contributions to the OPE of All exist 
already in the literature [8,18,19,24,38-40]. In some cases the quoted results (especially for 
as) differ significantly from ours. To pin down the source of these discrepancies, a closer 
scrutiny of the previous analyses is in order. In general, previous results have errors much 
larger than those on the spectral function over most of its measured range. This suggests 
either the impact of strong cancellations or the presence of additional theoretical systematic 
uncertainties. One obvious possibility is the presence of d > 8 contributions, neglected in 
the analyses of Refs. [8,18,19,24], in the solutions for a 6 , a 8 . We will demonstrate below 
that, for both and ag, the differences between our results and those of previous analyses 
are naturally accounted for by the d > 8 coefficients given in Eqs. (26), (28). 

In what follows, we shall recall the basic ingredients of the earlier analyses and discuss 
possible sources of uncertainty. 

A. Spectral Weight Analyses 

In Refs. [8,18,19], the "(k,m) spectral weights", 

w^ m \y) = y m (l-y) 2+k (l + 2y) , (29) 

with (k,m) = (0,0) and (l,m), m = 0, • • • , 3, were employed to extract a 6 and a§, under 
the implicit assumption that contributions with d > 8 were negligible in all cases 12 . The fits 
for a 6 and as were, in all cases, performed using only the highest So available, so = m 2 . 

Ref. [8] (DGHS) represents an update of the earlier ALEPH analysis [18], and con- 
centrates specifically on the V-A combination, which was not studied independently in the 
original ALEPH paper. The results for a$, a§ thus supercede those inferred from the separate 
V, A extractions performed in Ref. [18]. The results, in our notation, are 

a 6 = (-6.4 ± 1.8) x 10~ 3 GeV 6 

a 8 = ( 8.7 ± 2.4) x 10" 3 GeV 8 . (30) 
They are in good agreement with the results of the OPAL analysis [19], 

a 6 = (-6.0 ±0.1) x 10~ 3 GeV 6 

a s = ( 7.6 ± 0.6) x 10~ 3 GeV 8 . (31) 

One should bear in mind that the DGHS and OPAL analysis methods are somewhat differ- 
ent: the DGHS results follow from a dedicated V-A analysis, while the OPAL results were 
generated by combining the d — 6, 8 contributions extracted for the separate V and A corre- 
lators. The separate V, A analyses, however, involve an additional OPE fitting parameter, 



12 Ref. [8] also employed the (1, —1) spectral weight, not included in the other analyses, in order 
to allow the simultaneous extraction of the NLO chiral LEC L 10 . 
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the gluon condensate, which is absent in the V-A difference. The fits display very strong 
correlations between a^, a§ and the gluon condensate [19]. A dedicated V-A analysis of the 
OPAL data would thus, in general, be expected to give different results for a 6 , a 8 13 . In view 
of this, and the good agreement between the OPAL and DGHS results, we concentrate on 
the DGHS solution in the discussion which follows 14 . 

The DGHS value for a 6 is consistent with ours, within errors, but that for a$ is not. We 
have studied the origin of this discrepancy, and we find that: 

• The discrepancy can be understood as arising from the neglect of the d > 8 contribu- 
tions to the spectral weight sum rules employed by DGHS. 

• The "(k,m) spectral weights" FESR actually provide a consistency check on our solu- 
tion set. 

We first note that the (0, 0) and (1, 0) pFESR's have strong V-A cancellations, and hence 
large experimental errors on the data sides of the sum rules. For the (0, 0) case, which 
involves, from among the unknown d > 4 terms, only the and a 8 contributions, 
ry-A ~ 3% for So = m 2 .. The (1,0) case, whose OPE side in principle involves a^, a§ 
and aio, also has r V -A ~ 3% for s = m 2 .. The (1,1), (1,2) and (1,3) weights produce 
much less pronounced V-A cancellations 15 , and hence must dominate the DGHS fit. Note, 
however, that w^ 1,m " > (y) = y m [1 — y — 3y 2 + 5y 3 — 2y A ] ; these weights thus produce numerical 
enhancements of the a^, d > 8 terms, whose presence on the OPE sides of the sum rules 
has been assumed to be numerically negligible (see explicit example below). The pFESR's 
dominating the fit are thus those for which neglect of the d > 8 terms is least safe. 

It is easy to check that the combined fit values for aio, • • • , ai6 predict non-negligible d > 8 
contributions for all the (l,m) pFESR's. Our results thus imply that the DGHS values 
for ciq and a$, which are dominated by m = 1,2,3 cases, must contain higher dimension 
contamination. That the central DGHS a 6 , a 8 values do not provide as good a fit to the 



We thank Sven Menke for bringing this point to our attention. No analogue of the DGHS 
update of the ALEPH analysis exists, at present, for the OPAL data. 

14 A slightly different set of values, corresponding to an average of the results of Refs. [18] and 
[19], has been used in the (0,0) spectral weight analysis of Section 7 of Ref. [41]. The value for 
a 6 is the same as that of DGHS, while that for a§ is ~ 15% higher. The reader interested in the 
chiral limit value of the K — > -kit matrix element of the electroweak penguin operator, Q 8 , should 
bear in mind, not only the difference between the ag values of Refs. [8,41] and our results above, 
but also the fact that the extractions of the dominant, (Og), contribution to olq in Refs. [8,41] 
employ a value for the coefficient A% much larger than that given above. To convert {Og) as 
determined in Refs. [8,41] to the same renormalization scheme as used for the Wilson coefficients 
of the effective weak Hamiltonian (and hence to make meaningful comparisons with the results of 
Refs. [1,15,38,39,42]), one must multiply these results by factors 1.15 and 1.27 for the NDR and 
HV 75 schemes, respectively. 

15 For example, for the (1, 1) pFESR, r v _ A = 32% for s = m 2 T . 
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w\ and W2 pFESR's (for which d > 8 contributions are absent) as does our combined fit 
is, presumably, a reflection of this contamination. Further evidence is provided by the W3 
through wio pFESR's. 

One can also explicitly demonstrate that the neglected d > 8 contributions are, indeed, 
important for the spectral weight pFESR's. This demonstration is most transparent for the 
(1,3) pFESR since, in this case, the OPE integral is: 

, ,r -, \ a 8 OlO . Q Oi2 a 14 ai 6 

/ w( i,3) ({a d }; s ) = --3-- r + 3-^ + 5-^ + 2-^. (32) 

If d > 8 OPE contributions are indeed negligible then, rescaling J w (i,3) (s ) by Sq should pro- 
duce a result, — o 8 , independent of s . 16 We plot SoJ tu (i,3)(so) f° r the ALEPH data in Fig. 5. 
The result is clearly far from constant with respect to Sq, unambiguously demonstrating the 
presence of non-negligible d > 8 contributions. The solid line shows, for comparison, the 
predictions corresponding to the combined fit solution of Eqs. (26). The good match shows 
that the d > 8 contributions produced by our solution naturally account for the discrepancy 
between the DGHS predictions and the experimental results. 

A similar situation holds for the other spectral weight sum rules, as shown in Fig. 6. 
In the figure we display the J w (k,m)(so) together with the OPE expressions f w (k,m) {{ad}', so) 
corresponding to (i) the central values of the DGHS fit, Eqs. (30), together with a<j = for 
d > 8, (shown by the dashed line), and (ii) our combined (ALEPH-based) fit, Eqs. (26), 
(shown by the solid line). In all cases, if one takes into account the errors and correlation 
for the DGHS fit parameters, the resulting OPE error bar overlaps the spectral integral bar 
at s = m 2 ., even when the central values are not in particularly good agreement. However, 
when one goes to lower s this is no longer the case; the shape of the curve for the OPE 
integrals as a function of s is typically rather different from that for the spectral integrals. 
This is another signal of missing higher dimension contributions on the OPE sides of the 
sum rules. On the other hand when one considers the OPE contributions implied by our 
combined fit, a high-quality match between the OPE and data integrals is obtained. This 
is a non-trivial consistency test on our solution set. 



B. The IZ pFESR and Borel Sum Rule Analyses 

In Ref. [24] (IZ), three approaches were considered: (1) pFESR's with w(y) = (1 — y) 2 
and 2/(1 — y) 2 , (2) Borel transformed dispersion relations involving AU(Q 2 ) for Q 2 lying 
along various fixed rays in the complex Q 2 -plane, and (3) Gaussian sum rules. The results, 
in this case, are 

a 6 = (-6.8 ±2.1) x 1(T 3 GeV 6 

a 8 = ( 7 ± 4) x 1(T 3 GeV 8 , (33) 



16 This is valid up to small logarithmic corrections. For the DGHS solution, the correction as- 
sociated with the d = 6 logarithmic term varies from 2.3% to 1.5% as sq increases from 1.95 to 
3.15 GeV 2 . 
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and are dominated by the Borel sum rule (BSR) part of the analysis, though the other 
determinations are compatible with these, within their (larger) errors. 

The pFESR part of the IZ analysis involves one weight, w(y) = (1 — y) 2 , for which the 
V-A cancellation is extremely strong (r V -A ~ —4% for s = m 2 ), and one, w(y) = y(l—y) 2 , 
for which it is considerably less so (r V -A ~ 25% for s = rn 2 T ). The strong cancellation for 
the (1 — y) 2 case leads to large errors on clq, and to the strong sensitivity to the errors on 
F w noted in Ref. [24]. The necessity of subtracting the poorly determined contribution 
to the y(l — y) 2 sum rule before obtaining the residual a 8 contribution, then leads to large 
errors on a 8 as well. 

BSR's were employed in the second part of the analysis because of factorial suppression of 
high d contributions (a^ contributions appear in the Borel transform of the OPE side of the 
sum rule multiplied by l/((<i — 2)/2)\M d [43], where M is the Borel mass). Since, however, 
the spectral data are known only up to s = m 2 , and have significant errors above 2 GeV 2 , IZ 
are forced to work at quite low Borel masses to suppress contributions from the region of the 
spectrum where either data errors are large or data are absent. Explicitly, M 2 ~ 0.8 GeV 2 
is used for sum rules dominating the determination of a^, and M 2 ~ 0.6 GeV 2 in sum rules 
used for a$. At such low M 2 , factorial suppression of high d contributions is counteracted 
by the enhancement associated with the smallness of the M 2N+2 factor in the denominator, 
making the sum rules potentially sensitive to higher dimension contributions. 

While the central IZ values for a 6 and a 8 are obtained neglecting d > 8 contributions, 
the quoted errors include, not only the uncertainties due to experimental errors, but also a 
contribution meant to represent a plausible bound on the magnitude of the d > 8 terms. 
This bound is based on the assumptions that |aio| and | a.12 1 are bounded by 2 GeV 4 |ae| 
and 5 GeV 6 |a 6 |, respectively. According to the results of our fit, these assumptions are not 
sufficiently conservative: the bounds, in both cases, lie well outside the range allowed by 
the errors on the combined fit values. Recall also that, as shown in Figures 2, 3 and 4, 
the central IZ a^, a 8 values do not provide good fits to the w\ through wi$ pFESR's. In 
contrast, our combined fit implies values for the OPE sums for the four IZ BSR's which 
are in excellent agreement with experiment. A demonstration of this claim, together with a 
more detailed discussion of the four IZ BSR's, may be found in the Appendix. 

Comments similar to those on the BSR's apply to the Gaussian sum rules studied by IZ. 
Since, however, the Gaussian sum rule and as errors are larger than those of the BSR 
analysis, and the OPE convergence even slower, we will not comment further on that part 
of the IZ analysis. 



C. Duality Point Analyses 

A recent discussion of duality point analyses (summarized earlier in Section I) can be 
found in Ref. [9]. We comment here on the most recent numerical results, obtained in 
Refs. [38,39] (BGP) 17 . 



17 An estimate of the four-quark vacuum matrix elements which determine clq, obtained by trun- 
cating the spectral integrals appearing in the dispersive sum rules of Ref. [13] at the duality points 
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BGP determine ae and ag from FESR's based on the weights w(s) = s 2 and w(s) = s 3 , 18 
working at the highest duality point determined through the second WSR. From the analysis 
based on ALEPH data (for which = 2.53i°;i2 GeV 2 ) the following results are quoted: 

a 6 = - (3.4l|*) x 10- 3 GeV 6 

a 8 = - (l4.4^° 4 ) x 10~ 3 GeV 8 . (34) 

These values are in qualitative agreement with ours, but are affected by large uncertainties. 
The origin of these uncertainies is twofold. On the one hand the analysis uses weights which 
emphasize the region where the data errors are large and, on the other, the uncertainty 
in the exact location of the WSR duality point gets amplified by the strong slope of the 
relevant spectral integrals with respect to s near s = Sq. 

The errors on the second WSR duality point, quoted above are entirely experimental 
in origin. The fact that the duality points for the s 2 , s 3 FESR's may not coincide exactly with 
those of the second WSR, however, leads to an additional uncertainty which is not amenable 
to experimental improvement. We think this uncertainty is unlikely to be negligible, as 
argued below. If the duality points for the WSR's are universal, i.e. all FESR 's are satisfied 
at such Sq \ then the values of the extracted OPE parameters should not depend on the 
particular duality point used in the analysis. Empirically, however, if one uses the lower of 
the two second WSR duality points (s = 1.47±0.02 GeV 2 [38]), as advocated in Ref. [9,40], 
one obtains [38]: 

a 6 = - (13.2 ± 0.4) x 10~ 3 GeV 6 

a 8 = (24l 2 ) x 10- 3 GeV 8 . (35) 

These results have much smaller errors (reflecting the better data quality at lower s ) but 
are not compatible with those obtained using the higher duality point, Eq. (34). It is thus 
impossible for the two duality points of the second WSR to both be duality points of the 
s 2 and s 3 FESR's. Since at least one of the two WSR duality points must differ from the 
corresponding s 2 , s 3 duality point, it seems unlikely to us that either is exactly identical to 
its s 2 , s 3 counterpart. 



of the WSR's, was also given in Section 6 of Ref. [41]. The assumptions underlying this analysis 
are even stronger than those underlying the duality point truncation of the WSR's, where the 
corrections for the truncation can be shown to be numerically small. In addition, the original 
dispersive sum rule for the dominant (0$) contribution suffers from potential contamination by 
higher dimension effects [14] at the scale fi = 2 GeV employed in Ref. [41]. Since, in any case, the 
errors from this approach are a factor of > 3 larger than those obtained by averaging the results 
of the ALEPH and OPAL (0, 0) spectral weight analyses in Section 7 of the same reference, we do 
not discuss this estimate any further. 

18 A second determination of the dominant (Os) contribution to ae, which however requires an 
additional input assumption, gives a compatible result. 
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We emphasize that it is the strong slope of the data integrals with respect to sq which 
is particularly problematic for the duality point approach. The possibility of the existence 
of a reasonably narrow s region within which the actual duality points of a number of 
differently-weighted FESR's might lie is not itself implausible. Indeed, at those intermediate 
scales suggested by the PQW argument [12] (where OPE violation is small, except near the 
timelike real axis), duality points for a wide range of FESR's would be expected to cluster in 
the vicinity of any so for which the real and imaginary parts of n(so) — ^-ope(sq) happened 
to be simultaneously small. In the case of All, the zeros of Im [n(s ) — ^-ope(so)} on the 
real axis occur at s ~ 0.9 and 2.1 GeV 2 , somewhat removed from the locations of the 
WSR duality points. We would thus expect the s 2 and s 3 duality points to, indeed, differ 
somewhat from the corresponding WSR duality points. Since \Im [n(s ) — ^ope(so)] | is 
considerably smaller at the higher of the two WSR duality points, we are in agreement with 
the authors of Refs. [38,39] in expecting the higher of the two duality points to provide the 
more reliable estimate of a 6 and a 8 . This expectation would appear to be borne out by 
comparison to our results. In particular, the result that a 6 and a 8 have the same sign, first 
obtained in Ref. [38], is confirmed by our analysis. We remind the reader that the opposite 
sign for a$ obtained in both the spectral weight and BSR analyses is naturally accounted 
for by the d > 8 contributions implied by our solution but neglected in those analyses. 

D. The MHA Analysis 

In Ref. [40], the a^ >4 are determined, not from data, but using a large- iV c -inspired, 3- 
pole, model approximation to Ap (the so-called "minimal hadronic ansatz" or MHA). The 
s fe -weighted physical and MHA spectral integrals, for a given sq, are in general very different. 
For all —2 < k < 4, however, the point where the two agree happens to lie in the vicinity 
of the lower of the duality points for the two WSR's. This observation is taken as evidence 
in support of the pattern of long-/short-distance duality predicted by the MHA, and of the 
reliability of the model values of a 6 , a 8 and aio- 19 The results quoted in Ref. [40] correspond 
to 

a 6 = - (9.5 ± 2.0) x 10~ 3 GeV 6 
a 8 = (16.0 ± 4.2) x 10" 3 GeV 8 

a 10 = - (20.8 ± 10.2) x 10" 3 GeV 10 , (36) 

which are not in good agreement with our central fit values. 

One should bear in mind that the errors in Eqs. (36) reflect only the uncertainties in the 
fitted values of the three independent MHA parameters F , my, and 5% and not any possible 
theoretical systematic errors (due to working with the minimal set of hadronic states, and in 
the large N c limit). The latter are not necessarily negligible. As a first indication of this, let 
us observe that MHA predicts the duality points for the various s fc -weighted FESR's (those 



Recall that a,2k+2 = (— l) fc Jo* ^ s sk ^P( s ) ^ s o is a t rue duality point. 



20 



So for which the model and data integrals match) to be different 20 . Because of the strong 
slope of the s k (k = 2,3,4) spectral integrals with respect to Sq, even small errors in the 
model predictions for these differences can correspond to large uncertainties on the 

To further test the MHA predictions, and to get an idea of whether or not potential 
systematic uncertainties might account for the discepancy between the MHA predictions 
and our results, we may study those pFESR's sensitive only to a 6 , as and a w . We display, 
in Fig. 7, the MHA predictions for the J w (sq) associated with the w±, (0,0), and (1,0) 
pFESR's. (The first and second sum rules are sensitive to the MHA values of a^, a§, the third 
to a 6 and a 1Q , and the fourth to a 6 , a 8 , and a 10 .) We see that, although the representation of 
the physical spectral integrals is not-unreasonable for a three-parameter model, the quality 
of this representation is not good at the detailed level. This mismatch suggests to us the 
presence of residual theoretical systematics uncertainties in the MHA approach, which might 
be removed by going beyond the minimal ansatz and/or incorporating 1/N C corrections. 

V. DUALITY VIOLATION AND V-A PFESR'S 

Our interpretation of the ad obtained above as the true asymptotic OPE coefficients of 
the ud V-A correlator rests on the assumption that residual duality violation (parameterized 
by R[sq,w]) is small for the doubly-pinched pFESR's and scales employed in our analysis. 
As noted above, the high quality of the match between the optimized OPE representation 
and the corresponding spectral integrals is a necessary, but not sufficient, condition for the 
validity of this assumption. In this section we discuss additional evidence in its favor. We 
begin by reviewing certain relevant aspects of what is known about the nature of duality 
violation in QCD. 

A. General Expectations for Duality Violation in QCD 

A useful review of the current status of our understanding of duality violation in QCD 
is given in Ref. [17]. It is important to bear in mind that duality violation may be small in 
weighted spectral integrals even when the level of duality violation in the spectral function 
itself is large over significant portions, or even all, of the integration range 21 . 



20 This can be seen, for example, by superimposing the plots for the s° and s 2 moments (the top 
panels of Figs. 1 and 2 of Ref. [40], respectively). One finds that the band within which the s° 
matching point must lie (given the uncertainties in the model parameters) does not overlap with 
the corresponding allowed band for the s 2 moment. The situation is similar for the s 3 and s 4 
moments (the allowed matching bands in fact lies somewhat farther from the corresponding s° 
band than in the s 2 case). 

21 Examples are the spectral integrals corresponding to (i) dispersion representations of correlators 
for spacelike Q 2 3> Aq CD , (ii) Borel transformed dispersion relations involving Borel masses M » 
Aq CD and (iii) the (l—y) k (l+ Ay) -weighted (k = 1, 2) pFESR's for the flavor ud V and A correlators 
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Two distinct types of duality violation in the spectral function are identified in Ref. [17]. 
The first is that produced by contributions to the correlator which, asymptotically, are 
exponentially suppressed relative to OPE contributions for spacelike Q 2 . Such terms behave 
asymptotically as exp(— bQ)/Q K , and hence acquire oscillating imaginary parts for timelike 
Q 2 = —s. The second type of duality violation occurs in N c — > oo QCD, where the spectrum 
consists of a tower of infinitely narrow resonances. As a result of this spectral structure, 
the associated correlator has a different convergent Laurent expansion in each of the annuli 
lying between successive poles in the complex s-plane. In none of these annuli is the Laurent 
expansion equal to the asymptotic expansion; hence duality violation exists, in this case, in 
all such annuli. 

An important difference in the nature of duality violation in these two cases lies in the 
structure of the duality violating contributions to the correlator in the complex plane. In the 
N c — > oo scenario, one has a series of different "sub-asymptotic" expansions, each valid in a 
different annulus. When one crosses from one annulus to the next, all Laurent coefficients are 
altered, and for no annulus are they equal to the corresponding asymptotic OPE coefficients. 
Duality violation in a given annulus is thus equally large at all points on the circle \s\ = sq 
lying within the given annulus, and is NOT localized to the vicinity of the timelike point on 
that circle for any s , no matter how large. In contrast, for Q 2 = s e 1 ^ (with = — n (+7r) 
corresponding to the top (bottom) of the physical cut), a term of the form ~ exp(— bQ)/Q K 
behaves as 

\ exp [-b s cos(0/2)] exp [i («0/2 + sin(0/2))] (37) 

and hence retains an exponential suppression, via the factor exp[— sob cos(0/2)], for all but 
the timelike point on \s\ = s . This suppression will remain quite significant over most 
of the circle for scales s larger than 1/b. This may be the case even if the oscillating, 
duality- violating component of the spectral function is far from negligible at the same sq. 
Duality violation via such terms is thus PQW-like: "intermediate" scales exist for which 
duality violation in the correlator (IIope(s) ~ n(s)) is strongly localized to the vicinity of 
the timelike real axis. 

B. Detecting the Presence and Nature of Duality Violation 

The two distinct patterns of duality violation for a given correlator in the complex 
plane manifest themselves in readily distinguishable different ways in sum rule analyses. 
In particular, these patterns suggest different strategies and tests to explore the impact 
of duality violation in a given analysis. In this section we identify such strategies and 
enumerate a number of possible tests. In the next section, we then specialize to the flavor 
ud V-A correlator, and we discuss the practical implementation of these tests. 

In presence of a PQW-like component of duality violation, there will exist intermediate 
scales s where 



at scales so ~ 2 to 3 GeV 2 . 
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• power-weighted (w(s) = s k ) FESR's, which fail to suppress contributions from the 
integral over \s\ = s near the timelike axis, are poorly satisfied; 

• pFESR's involving weights which suppress contributions from the vicinity of the time- 
like point are well satisfied. 

This observation motivates the use of "pinched weights" in FESR analyses to tame PQW- 
like duality violation. Having adopted a set of weights, one wants to verify that residual 
duality violating contributions from the region near s = Sq are not present in the results of 
the analysis. In this respect, an important test is to 

1. verify that the (nominally asymptotic) coefficients extracted in the pFESR analysis 
provide accurate representations, not only of the spectral integrals used in fitting the 
but also of spectral integrals corresponding to weights with zeros of a higher order 
at s = so (which therefore further suppress PQW-like duality violating effects). 

In the N c — > oo-like scenario, where duality violation is not localized to the vicinity 
of the timelike real axis, a rather different pattern of sum rule behavior will be observed. 
Since OPE-like Laurent expansions exist in any given annulus, so long as one restricts 
oneself to Sq lying in a single annulus, one will obtain a set of coefficients, a^, which provide 
a perfect match between the "OPE" and data sides of both power-weighted and pinch- 
weighted FESR's at those scales. That set will, of course, consist of just the coefficients of 
those terms in the Laurent expansion for the given annulus which survive when integrated 
against the weights employed. If one performs the same basic analysis (i.e., using the same 
set of weights), but now for s lying entirely in a different annulus, one will obtain a different 
set of ad- These ad will provide a perfect match between the "OPE" and data sides of the sum 
rules employed in the new annulus. Schematically, this type of duality violating contribution 
implies: 

• the existence of several sub-asymptotic regimes; 

• that pinching is not effective in removing this type of duality violating effect (both 
pinch-weighted and power-weighted FESR's are equally well satisfied in each sub- 
asymptotic region). 

The existence of several such sub-asymptotic regimes in this type of scenario can also, in 
principle, be exposed in a pFESR analysis, as follows. Starting with some particular small 
range of s one ma Y gradually decrease the lower edge of the s analysis window, keeping 
the upper edge fixed. So long as the lower edge lies in the same annulus as the upper edge, 
a pFESR analysis extraction of OPE-like coefficients, obtained by means of matching to 
spectral integral data, will produce an exact determination of the relevant coefficients in 
the singular part of the Laurent expansion for that annulus. As soon as the lower edge 
of the s analysis window leaves the single annulus, however, there is no longer a single 
set of expansion coefficients valid for all the s being employed. The existence of such a 
two-annulus regime will be evident by the sudden appearance of a poor match between the 
optimized OPE-like integral and spectral integral sets. One can check these expectations 
explicitly within the equal-spacing pole model of Refs. [9,17]. This exercise shows that, in 
performing a pFESR analysis, it is crucial to 
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2. demonstrate that there is no drift in the values of the extracted coefficients as one 
decreases the lower edge of the s analysis window and 

3. demonstrate that the optimized values of the fit parameters in fact produce an accurate 
match to the spectral integral data used to produce the fit over the whole of the s 
analysis window employed. 

These tests serve not only to verify the reliability of the assumed OPE-like expansion form, 
but also, at least potentially, to expose the existence of multiple sub-asymptotic expansion 
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regimes . 

In the context of the N c — > oo discussion, however, it is clear that, while passing these 
tests is a necessary condition for the reliability of the extraction of the OPE coefficients 
from data, it is not a sufficient one: if one happened to be unlucky and perform the pFESR 
analysis only for those So lying in a single, but sub-asymptotic, annulus, one would see a 
high quality match (exact in the case of the pole model) between the spectral integrals and 
optimized OPE-like integrals even though one would have actually extracted the coefficients 
relevant to the Laurent expansion in the sub-asymptotic annulus, and not those relevant to 
the asymptotic regime. A simple way to test whether or not this is the case is to 

4. take the coefficients extracted in the pFESR analysis and employ them as input to a 
dispersive analysis relevant to the asymptotic regime. 

If the coefficients extracted in the pFESR analysis are not those relevant to the asymptotic 
regime, the resulting dispersive integrals will be poorly approximated by the OPE-like repre- 
sentation generated using the fitted coefficients. Again, explicit illustrations of this point can 
be worked out within the equal-spacing pole model of Refs. [9,17]. In the case of the model, 
performing the dispersive test is straightforward because the spectral function of the model 
is actually known for all s. The situation of interest to us, however, is one where spectral 
data are available for only a limited range of s. In such a situation, it would typically be diffi- 
cult to construct a dispersive test for which the errors on the dispersive integrals were under 
sufficient control to make the test useful. One general solution is to work with BSR's and 
restrict one's attention to Borel masses which are both low enough that the spectral weight, 
exp(— s/M 2 ), is negligible in the region where spectral data are absent and, simultaneously, 
high enough that the convergence with d of the Borel transformed OPE series for the set of 
ad one wishes to test is acceptable. For a given BSR, such M may or may not exist. Three 
of the four IZ BSR's turn out to provide examples of such tests for our solution set (details 



One may also test whether an observed deterioration in the quality of the optimized 
"OPE" /spectral integral match as sq is lowered is, or is not, due to the existence of a new sub- 
asymptotic regime. If it is, then the first so for which the deterioration appears must lie in the 
lower annulus. Working with sq lying in a narrow range just below this point should then produce 
a new set of fitted ad which provide a good quality representation of the corresponding spectral 
integrals, when restricted to this new range of sq. If a good quality match is not found, then the 
deterioration is due to breakdown of the OPE-like expansion form, and not to the fact that one 
has entered a new subasymptotic region. 
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are reported in the Appendix). Additional asymptotic tests, involving BSR's at larger M, 
are possible for the flavor ud V-A correlator. In this case, the spectral integral uncertainties 
are brought under control using the classical chiral sum rule constraints associated with the 
Weinberg sum rules and the sum rule for the ir electromagnetic mass splitting. An efficient 
procedure for implementing these constraints is provided by the "residual weight method" , 
which is described in detail in Ref. [15]. 

C. The Nature of Duality Violation in the ud V-A Channel 

In the following we argue that duality violation in the ud V-A correlator is predominantly 
PQW-like. From previous work [16], one has empirical evidence that duality violation in 
the individual ud V and A channels is predominantly PQW-like. Checking for the presence 
or absence of duality violations in the ud V and A correlators at intermediate scales is 
straightforward because one has independent (asymptotic) information on the value of the 
OPE parameter a s . Such a straightforward check is not possible for the ud V-A difference. A 
qualitative argument is however available, based on the the observation that duality violation 
in the flavor ud V+A sum cancels, within experimental errors, for scales above s ~ 2 GeV 2 . 
This can be seen from (i) the fact that the corresponding spectral function is in agreement 
with the OPE prediction for such s (see, e.g., Fig. 6 of the second of Refs. [18]) and (ii) the 
observation that the spectral integrals for the s fe -weighted FESR's, k — 0, • • • , 3, are in good 
agreement with the corresponding OPE integrals for sq above ~ 1.9 GeV 2 [16]. This implies 
that the duality violating contribution to the ud V-A correlator is, within experimental 
errors, twice that to the ud V correlator. The latter is known to be strongly localized to the 
vicinity of the timelike real axis for the scales of interest to us, and hence so is the former. 
This conclusion is compatible with the observation that, although the ud V-A s fc -weighted 
FESR data integrals are not constant with respect to so (i.e., not in agreement with the 
behavior of the s fc -weighted "OPE" integrals), the agreement between the data and "OPE" 
sides of our pFESR's is very good for the optimized OPE-like fits given above. 

Tests of the type 1. 

Having a suppression of duality violating contributions which is strong enough to make 
such contributions negligible relative to the d = terms in the ud V and A correlators 
does not necessarily mean that the same suppression is sufficient to make such contributions 
small relative to the d = 6 and higher OPE contributions in the ud V-A difference. In order 
to check for residual duality violating contributions localized to the vicinity of the timelike 
real axis, we have performed tests of the type defined in the previous section, item 1. The 
(1, to) spectral weights discussed above (with to = 0, • • • , 3) have zeros of order 3. As we 
have already seen in Figs. 5 and 6, the results of our combined fit produce an extremely 
good "OPE" /spectral integral match for all of these weights, with no quality deterioration. 
We have also investigated the (2,0), (2, 1), (3,0), (3, 1) and (4,0) spectral weight pFESR's, 
which have weights with zeros of order 4, 4, 5, 5 and 6, respectively, at s = s . Again 
the quality of the match to the spectral integral sides of these sum rules provided by our 
combined fit is excellent in all cases, despite the much stronger suppression of contributions 
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from the region on \s\ = Sq near s = Sq. We illustrate the quality of this match for the most 
extreme cases (the (3, 1) and (4,0) pFESR's) in Fig. 8. 

Tests of the type 2. ,3. ,4. 

The arguments given above do not completely rule out the presence of residual duality 
violation of non-PQW-like nature (the N c — > oo-like scenario). In order to deal with this, 
we have subjected our solution set to tests of the type described in items 2. ,3., and 4. of the 
previous section. As for test 2., we find that within the present experimental errors there 
is no drift in the extracted OPE parameters as one lowers the lower edge of the So analysis 
window (see below for details and prospect of sharpening this test with improved data). 
Also tests of the type 3. are succesfully passed by our solution set (see Section III C) . 

Finally, to deal with the possibility that our entire s analysis window lies within a single 
sub-asymptotic region, we have performed a number of asymptotic dispersive tests of the 
type described in the previous subsection, item 4. A first set of asymptotic dispersive tests is 
provided by the four IZ BSR's. These are highly non-trivial since, because of the difference 
in the sign of ag between our combined fit and the IZ solutions, those IZ BSR's for which 
d — 6 contributions are absent would appear to be problematic for our combined fit. It turns 
out that this is not the case; in fact, the convergence of the Borel transformed OPE series 
is quite slow at the low M employed by IZ and, once one extends the sum involving our 
combined fit to sufficiently high d to obtain convergence, the OPE predictions are in excellent 
agreement with the spectral data. Since these tests are also relevant to the comparison to 
previous work, we provide a detailed demonstration of these claims in the Appendix. 

To obtain BSR's at larger Borel mass, M, one needs to use "residual weight method" 
improvement on the spectral integrals [15]. In order to keep the errors under control, it 
is necessary to work with the product of All with appropriately chosen polynomials. We 
find that the combined fit OPE predictions are in excellent agreement with the spectral 
integral sides of these BSR's for M over a range sufficiently wide that, at the upper end, 
the OPE integrals are completely dominated by their d — 6 contribution while, at the lower 
end, the full set of obtained in the combined fit (d — 6, • • • , 16) must be included before 
convergence of the Borel transformed OPE sum is obtained. 

Model explorations 

In principle, explicit models of the V-A spectral function could be used to try and address 
the level of duality violation present in our analysis. One should bear in mind, however, 
that the only information we have about Ap in the region above s = m^. is in the form of 
the constraints provided by the classical chiral sum rules. These constraints are far from 
sufficient to fully constrain the behavior of Ap above s = and, as a result, there exists 
a wide range of model extensions of the data for Ap to s > m%., all of which are compatible 
with these constraints. It is easy to construct, among these, models for Ap(s) for which the 
asymptotic expansion coefficients are the same as those of our combined fit. The models 
which have this property display continued damped oscillations in Ap as one goes to higher 
s, and hence appear quite natural. It is also possible to construct models for which the 
asymptotic OPE parameters differ significantly from those of our combined fit [44]. 



26 



Because the integrated pFESR OPE contributions of dimension d scale as 1/sq , one 
finds that, for large Sq, the higher d ad contributions drop rapidly in size with increasing d. 
With such small high d contributions, a small change in the modelling of Ap in the region 
where it is not known experimentally typically produces a large change in a^- A very large 
theoretical systematic uncertainty for the higher dimension will thus be associated with 
any attempts to model Ap in the region above s = m 2 . Without being able to control 
this theoretical systematic error, obtaining meaningful information on the level of duality 
violation from such model studies is somewhat problematic. 

Prospects of improving the data-based tests 

It is worth stressing that significant improvements in the analysis will become possible 
once the new hadronic r decay data from the B-factory experiments is available. At present 
both the errors on the and the accuracy with which it is possible to determine the location 
of the onset of duality violation in the analysis are limited by the errors on Ap(s) above 
s ~ 2 GeV 2 . These errors are dominated by experimental uncertainties on the 4n, KKn and 
KKnn spectral distributions and uncertainties in the V/A separation for KK7r and KK7nr 
states. Major improvements should be forthcoming as a result of the expected ~ 10 2 -fold 
increase in the size of the r decay data base. The improved spectral integral errors which 
result will allow us to improve significantly on the efficiency of our tests for the absence of 
residual duality violation. The current situation in this regard is discussed in brief below. 

Recall that, by decreasing the lower edge of the analysis window, we were able to demon- 
strate the presence of duality violation for the pFESR's used in our analysis at scales below 
~ 1.8 GeV 2 . With current experimental errors there is no evidence for duality violation in 
our analysis window. Ideally one would like to work at scales well above 1.8 GeV 2 , in order 
to suppress, as much as possible, any residual duality violating contributions which might be 
present, but masked by current experimental errors. While current errors are small enough 
that a 6 may still be determined, even if one works with only a small portion of our present 
analysis window 23 , this is not true for the aa with d > 8. In fact, with current experimental 
errors, the uncertainties on the extracted d > 8 do not become smaller than |a<z|, until 
the lower edge of the analysis window has been reduced to below so ~ 2.5 GeV 2 . If, as an 
example, we perform our analysis of the ALEPH data using the Sq = 2.35 — > 3.15 GeV 2 
sub- window, then, with central values for all non-spectral input, the maximally-safe output 
for a 6 and a 8 is: 

a 6 = -(3.88 ± 1.21) x 1(T 3 GeV 6 

a 8 = -(9.32 ± 6.58) x 10~ 3 GeV 8 . (38) 



23 For example, using only sq = 2.75,2.95 and 3.15 GeV 2 , one finds, from the maximally-safe 
analysis of the ALEPH data, = —0.0049 ± 0.0022, where the error quoted is that associated 
with the ALEPH covariance matrix. The error is, of course, significantly larger than that obtained 
from the larger analysis window, but still less than 50% of the signal. 
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Within the quoted errors, these results are compatible with those of the full-window analysis. 
The situation is similar for the results of the combined analysis: the results of the sub- window 
analysis for a w through a i6 are: 

a w = (6.62 ±2.83) x 1(T 2 GeV 10 
a 12 = -(2.16 ± 0.86) x KT 1 GeV 12 
a 14 = (5.88 ± 2.48) x KT 1 GeV 14 

a 16 = -(1.47 ± 0.69) GeV 16 , (39) 

again compatible with the full-window analysis within the sub-window analysis errors. Were 
the errors to be 1/3 as large, however, the full- window and sub- window results would no 
longer be compatible and we would be forced to conclude that residual duality violation was 
present for those sq in the lower part of the full analysis window. We stress that there is 
no reason for reaching such a conclusion at present. In fact, there are strong reasons for 
trusting the results of the full-window analysis: 

• where the existence of duality violation can be explicitly demonstrated, the OPE-like 
expansion is known not to provide a good representation of the spectral integrals; 

• in the lower part of our full analysis window the OPE-like form provides an excellent 
representation of the spectral integrals; 

• the combined fit from the full analysis window provides an excellent representation 
of the spectral integrals not only in the lower part, but also the upper part, of the 
analysis window; 

• the combined fit results obtained from the sub-window version of the analysis turn out 
to provide a poor representation of the spectral integrals in the lower part of the full 
window. 

Nonetheless, the size of the sub-window errors are such that much stronger tests of the ab- 
sence of residual duality violation, using various sub-windows, will become possible once the 
errors on Ap(s) above 2 GeV 2 are reduced 24 . While it seems unlikely to us, for the reasons 
given above, it is not at present possible to conclusively rule out additional uncertainties, 
associated with residual duality violating, at the level of the difference of the full-window 
and sub-window analysis central ad values. 

VI. CONCLUSIONS 

In this paper we have used Finite Energy Sum Rules with "pinched-weights" (pFESR's) 
to determine the OPE coefficients ae,---,aiQ of the flavor ud V-A correlator with good 



If such reduced errors were to expose residual duality violation in the lower part of the current 
full- window analysis, one would of course be forced to raise the lower edge of the analysis window. 
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accuracy using existing hadronic r decay data. While it is not possible at present to either 
prove or disprove on rigorous analytic grounds that this approach (or any other) yields a 
valid approximation to the actual dynamics of QCD, we have carefully demonstrated the 
advantages of pFESR's among the class of sum rule techniques and have described a large 
number of checks on our own work and that of others. 

At a technical level, we have employed a set of ten polynomial weights carefully chosen 
to minimize the impact of experimental errors and of duality violating effects, as well as to 
optimally separate the contributions from condensate combinations of different dimension. 
Our analysis shows that the OPE contributions with d > 8 are typically not negligible at 
scales ~ 2 - 3 GeV 2 . 

We have performed a number of tests to explore the presence of duality violating effects 
in our analysis. These support the conclusion that our combined fit values are not affected 
by duality violation within the existing experimentally induced errors. We recall the main 
observations in support of this statement: 

(i) independent determinations of the using pFESR's based on different (independent) 
weights are in excellent agreement; 

(ii) the results of the combined fit for the ad lead to an extremely good match between the 
OPE and spectral integral sides of all the pFESR's employed in the fitting procedure; 

(iii) the combined fit values also lead to extremely good matches for the (k, m) spectral 
weight pFESR's, where d > 8 contributions are much larger relative to d = 6,8 
contributions than is the case for the W3 through wio pFESR's; 

(iv) there is no deterioration in the quality of the combined fit prediction for the pFESR 
spectral integrals even for those spectral weights with zeros at s = s °f much higher 
order than those used in obtaining the combined fit ; 

(v) the dispersive tests, described above, and in the Appendix, are successfully passed by 
our solution set. This provides additional support for the reliability of the extracted 
values, and our interpretation of them as asymptotic OPE coefficients of the V-A 
correlator. 

Improved experimental data would allow one to significantly sharpen some of the tests 
reported above. 

Some general observations also follow from the results and discussion above. First, the 
OPE representation of An, with the given by the combined fit values of either Eq. (26) 
or Eqs. (28), provides a very accurate representation of the corresponding spectral integrals 
down to scales as low as so — 2 GeV 2 , at least for pFESR's based on weights with a 
double zero at s = s . This suggests that the OPE remains reliable at intermediate scales, 
Q 2 ~ 2 — > 3 GeV 2 , apart perhaps from a region near the timelike real axis. In contrast, if 
one considers weights which do not suppress contributions from this region, one sees clear 
evidence for the breakdown of the OPE. The situation is similar to that for the flavor ud 
V and A correlators. The double zeros of the pFESR weights at s = sq in the V-A case 
evidently again provide sufficiently strong suppression in the vicinity of the timelike real axis 
to efficiently remove contributions from the region of OPE breakdown on the circle \s\ = Sq. 
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A second point concerns the relative sizes of the various ad- The results of the combined 
fit indicate that a^/^d is typically of order 2-^3 GeV 2 for the V-A correlator. This means 
that, at intermediate (2^3 GeV 2 ) scales, there is no 'natural' ordering of contributions 
with different d, in the sense that pFESR weights with comparable coefficients for the y N 
and y M terms in w(y) will produce comparable d = 2N + 2 and d = 2M + 2 contributions to 
the pFESR OPE integrals. This makes explicit the danger of neglecting terms with d > 8 
for pFESR's based on w(y) with degree greater than 3. This observation also raises the pos- 
sibility that the analogous neglect of higher d terms in other pFESR analyses, such as those 
used to extract m s from the flavor-breaking difference of ud and us V+A correlators [26,45], 
may suffer from similar problems. 
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APPENDIX A: THE IZ LOW-SCALE ASYMPTOTIC DISPERSIVE TESTS 

In this Appendix we provide details of the four IZ BSR's and complete the comparison of 
our results to those of early work by subjecting our combined fit to the asymptotic dispersive 
tests provided by these sum rules. 

Incorporating the small d = 6 logarithmic contribution, the BSR's employed by IZ may 
be cast into the form 25 



25 In writing Eqs. (Al), the factor F% from the RHS of Eq. (21) of IZ has been moved to the LHS 
of Eq. (Al), and absorbed into the spectral function via the shift Ap^ — ► Ap( 0+1 \ Eq. (A2) is 
just M 2 times Eq. (22) of IZ, up to the logarithmic correction term (proportional to bo). 
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(A2) 



where is the angle fixing the ray in the complex Q 2 plane along which the Borel transform 
was performed (0 = corresponds to the top of the physical cut). The four cases considered 
by IZ correspond to 

1. Eq. (Al) with = 5tt/6 and M 2 = 0.8 GeV 2 , 

2. Eq. (A2) with = 2tt/3 and M 2 = 0.85 GeV 2 , 

3. Eq. (Al) with = 3tt/4 and M 2 = 0.6 GeV 2 , 

4. Eq. (A2) with = 3tt/4 and M 2 = 0.65 GeV 2 . 

The first two cases have no d = 8 contribution, the third no d = 6 contribution. We test the 
combined fit by employing the fitted ad values as input on the OPE side of the IZ BSR's. 
This leads to a prediction for the value of the corresponding spectral integral (less the known 
d — 4 contribution) for each such sum rule. 

Since the a$ or a§ values obtained by IZ reflect the values of the spectral integrals, 
the change in sign of a 8 between the IZ fit and our combined fit would seem to represent 
a problem for the combined fit, especially in the case of the third IZ sum rule. This is, 
however, not the case. It is easy to check that, with the combined fit values as input, the 
convergence of the Borel transformed OPE series, at the low values of M 2 employed by IZ, 
is rather slow. Not only is the first of the d > 8 contributions neglected by IZ in obtaining 
their central values, in all cases, larger in magnitude than the corresponding sum of d — 6 
and/or d — 8 contributions, but also, in order to be certain that we have reached the region 
of convergence, we have had to extend the extraction of the ad to higher d. This is done 
using the extensions to higher N of the families of weights of which wq and wio are members 
(recall that these weights produce d > 4 OPE contributions proportional to either a 6 and 
02^+4 or as and 02.^+4) ■ We are able to extract terms with d up to 24, albeit with larger 
errors than for the given by the combined fit. The d = 18, •••24 values obtained 
from the two different weight families are in good agreement, and the a 6 and a 8 values 
obtained using each of the new sum rules separately are also in good agreement with those 
corresponding combined fit values. With these values in hand one finds that the OPE sides 
of the BSR's may be safely truncated, as can be seen from Table I. As in turns out, the 
additional (d > 16) terms play a role in the full OPE sum only for the third of the IZ sum 
rules. In none of the four cases is the d = 6 or d = 8 contribution the dominant one. 

In Table II we show the combined-fit predictions, together with the actual values (inferred 
from IZ) for the IZ spectral integrals (less the known d = 4 OPE terms). Shown for 
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comparison are the predictions corresponding to the IZ solution. Note that the second of 
the four sum rules has been used by IZ to obtain their quoted value for a^. The errors quoted 
for the combined fit d = 6 through d — 24 sum are obtained using the covariance matrix for 
the solution set, generated from covariance matrix of the ALEPH data 26 . As is evident from 
the table, the combined fit predictions are in excellent agreement with the data, providing 
further support for the asymptotic nature of the coefficients extracted in the combined fit. 
Table I, in addition, shows that the four sum rules weight the different ad contributions 
in very different ways, demonstrating that the results represent four independent, highly 
non-trivial tests of the combined fit results. 



IZ case 


d = 6 


d = 8 


d= 10 


d= 12 


d= 14 


d= 16 


d= 18 


d = 20 


d = 22 


d = 24 


1 


1 





1.33 


-1.93 


1.23 


-0.46 


0.09 


0.00 


-0.01 


0.00 


2 


1 





1.18 


-0.93 


0.00 


0.20 


-0.07 





0.00 


0.00 


3 





1 


-5.19 


4.09 


0.00 


-1.73 


1.16 


-0.34 


0.00 


0.04 


4 


1 


0.42 


0.00 


-1.46 


1.41 


-0.53 





0.09 


-0.04 


0.01 



TABLE I. The relative size of d > 4 OPE contributions to the four IZ BSR's for the extended 
version of the combined fit described in the Appendix. In all cases, the entries have been normalized 
to the lowest dimension (d = 6 or d = 8) contribution. IZ cases 1 through 4 label the four IZ BSR's 
according to the enumeration scheme given in the Appendix. 



"Because there are significant cancellations amongst the d > 4 contributions, and strong corre- 
lations among the combined fit values for the a^, it is crucial to employ the full covariance matrix 
for the solution in determining the uncertainties on the combined fit predictions shown in Table II. 
The fact that, at the low values of M for which the spectral integral errors are under control, the 
OPE sides of the IZ BSR's turn out to involve a complicated cancellation between a large number 
of terms of different dimension, in fact means that it would not be possible to disentangle such 
contributions using a BSR analysis. 
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IZ case 



Combined fit 



IZ fit 



Data - (d = 4) 



2 
3 
4 



-0.0023 ± 0.0002 
0.0038 ± 0.0002 

-0.0030 ± 0.0004 
0.0050 ± 0.0001 



-0.0027 ±0.0008 
0.0041 ±0.0013 

-0.0038 ± 0.0022 
0.0050 ± 0.0033 



-0.0023 ± 0.0006 
0.0041 ± 0.0009 

-0.0032 ± 0.0009 
0.0051 ± 0.0004 



TABLE II. The combined fit d > 4 predictions for the four IZ BSR's. IZ cases 1 through 4 
label the four IZ BSR's according to the enumeration scheme given in the Appendix. Column 2 
gives our prediction for the d > 4 OPE sum, obtained using extended version of the combined 
fit corresponding to the ALEPH data. The results obtained using the OPAL data are the same, 
except in the fourth case, where our prediction becomes .0051 ± .0001. Column 3 gives the d > 4 
OPE sum corresponding to the central values of the IZ fit. Column 4 gives the spectral integrals, 
less the known d = 4 terms, corresponding to the results quoted by IZ. The entries in columns 2, 
3, and 4 are in GeV 2 . 
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FIG. 1. The ALEPH (top panel) and OPAL (bottom panel) versions of the V-A spectral 
function. The errors shown are the square roots of the diagonal entries of the corresponding 
covariance matrices. 
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FIG. 2. J Wn (so) and fw„({ad}', so) for the w\ (top panel) and W2 (bottom panel) pFESR's. 
The J Wn (so) integrals and errors were obtained using the ALEPH data and covariance matrix. 
Three versions of the f Wn ({ad}', sq) curve are shown. The solid line corresponds to either the 
"maximally-safe" or combined fit for a§ and as, as described in the text, the short-dashed and 
long-dashed lines to the corresponding DGHS and IZ solutions, respectively. 
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FIG. 3. J Wn (so) and f Wn ({ad}; so) for the W3 (top left panel), w/4 (top right panel), w$ (bottom 
left panel) and wq (bottom right panel) pFESR's. The J w „(so) integrals and errors were obtained 
using the ALEPH data and covariance matrix. Three versions of the f Wn ({dd}', s o) curve are shown. 
The solid line corresponds to the combined fit given in Eqs. (26) of the text, the short-dashed and 
long-dashed lines to the DGHS and IZ solutions (for which aa = for d > 8). 
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FIG. 4. J Wn (so) and f Wn {{o-d\'i sq) for the wj (top left panel), w§ (top right panel), wg (bottom 
left panel) and w\o (bottom right panel) pFESR's. The J Wn (so) integrals and errors were obtained 
using the ALEPH data and covariance matrix. The notation for the three OPE curves is as in 
Figure 3. 
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FIG. 5. The rescaled (1,3) spectral weight combination, Sq J w {i,3) {sq) versus Sq. The integrals 
J w (i,3) (sq) and errors were obtained using the ALEPH data and covariance matrix. The solid line 
shows the OPE prediction f w (i,3)({ad}; sq) corresponding to the combined fit of Eqs. (26). 
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FIG. 6. J w (k,m)(so) and f w (k,m){{ad}; sq) for the (0,0) (top left panel), (1,0) (top right panel), 
(1, 1) (bottom left panel) and (1, 2) (bottom right panel) spectral weight pFESR's. The J w (k, m ) (so) 
integrals and errors were obtained using the ALEPH data and covariance matrix. The dashed and 
solid curves correspond to the OPE fit of DGHS, and our combined fit (Eqs. (26)), respectively. 
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FIG. 7. J w (so) and the corresponding MHA integrals for the w\ (top left panel), (0,0) (top 
right panel), W3 (bottom left panel) and (1,0) (bottom right panel) pFESR's. The J w (so) integrals 
and errors were obtained using the ALEPH data and covariance matrix. 
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FIG. 8. J w (k,m)(s ) and f w (k, m )({a d }; s ) for the (3,1) (top panel) and (4,0) (bottom panel) 
spectral weight pFESR's. The J w (k,m)(so) integrals and errors were obtained using the ALEPH 
data and covariance matrix. The solid line shows the predictions based on the d > 4 OPE solution 
represented by our combined fit, Eqs. (26). 
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